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ABSTRACT 

In the excursion set approach to structure formation initiahy spherical regions of the 
hnear density field collapse if the average density contrast within them exceeds some 
critical value, 5c- This model allows one to make predictions for several fundamental 
aspects of structure formation such as the mass function of dark matter halos, their 
merger histories and clustering properties. The value of 5c is often calculated from the 
spherical or ellipsoidal collapse model, which provide well-defined predictions given 
auxiliary properties of the linear tidal field at a given point. We use two cosmological 
simulations of structure growth in a LCDM cosmology to test a key assumption used 
in calculating 5c'- that the shapes of the initial Lagrangian patches that eventually 
collapse (or proto-haloes) are spherical. Our results indicate that the vast majority 
of dark matter proto-haloes are non-spherical, and have minor-to-major axis ratios 
that vary from (03/01) ~ 0.4 at the galaxy mass scale to ^ 0.65 for rich galaxy 
clusters. We show that this non-sphericity likely originates from the asymmetry of 
the linear tidal field which pushes material onto, or away from, local density maxima 
in the linear density field. We study the implications of these results for the collapse 
barriers for CDM halo formation inferred from the classic ellipsoidal collapse model. 
Our results indicate that the "standard" ellipsoidal collapse model commonly adopted 
in the literature does not provide a full account of the possible collapse thresholds for 
halo formation, since the model predictions depend sensitively on the assumed shape 
of the primordial perturbation. We show that an improved model, which accounts 
for the intrinsic shapes of proto-haloes, provides a much more accurate description of 
their measured minimum overdensities. 
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1 INTRODUCTION 

In standard models for the growth of large scale structure, 
such as the A-cold dark matter (CDM) model, dark mat- 
ter haloes form hierarchically from low-amplitude Gaussian 
density fluctuations seeded during inflation. Galaxy forma- 
tion is assumed to take place as gas cools and condenses 
in the centers of these haloes which subsequently "light-up" 
the underlying large scale matter distribution (e.g., White 
& Rees, 1978). Although the hydrodynamical processes in- 
herent to the formation of galaxies are difficult to model 
reliably, they are largely irrelevant to the over-all formation 
and clustering of the dark matter haloes. A flrst step to- 
wards understanding the statistics of the galaxy distribution 
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is thus a detailed understanding of the spatial distribution 
and clustering of their haloes. 

Since CDM haloes are highly non-linear objects much 
of our understanding of their properties has come from di- 
rect N-body simulations of their formation. Although use- 
ful, these simulations are often difficult to interpret and our 
understanding of the structure formation process can be sig- 
nificantly enhanced by augmenting them with simple phys- 
ical models for the collapse and virialization of dark matter 
haloes. 

One successful analytic approach to modeling struc- 
ture formation has become known as the Extended Press- 
Schechter formalism (see, e.g.. Press & Schechter, 1974; 
Bond et al., 1991). Here one identifies collapsed objects at 
a given time with regions of the linear density field expected 
to collapse based on some dynamical model, such as the 
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spherical (Gunn & Gott, 1972) or ellipsoidal collapse model 
(e.g., Zel'Dovich, 1970; Icke, 1973; White & Silk, 1979; Hoff- 
man, 1986; Lemson, 1993; Bertschinger & Jain, 1994; Bond 
& Myers, 1996). Although Press & Schechter (1974) orig- 
inally used this model to calculate the abundance of col- 
lapsed objects as a function of mass and redshift, it can 
be extended to compute a variety of hierarchical clustering 
statistics, such as halo merger histories (e.g., Lacey & Cole, 
1993; Sheth & Lemson, 1999b; van den Bosch, 2002), and 
spatial correlations (Mo & White, 1996; Mo et al., 1997; 
Catelan et al., 1998; Sheth & Lemson, 1999a). A related 
approach, the peaks formalism, relates the abundance of 
bound objects to the abundance of local density maxima 
above some threshold in the smoothed linear density field 
(Peacock & Heavens, 1985; Bardeen et al., 1986; Ludlow & 
Porciani, 2011). 

A fundamental aspect of both the Extended Press- 
Schechter approach and the peaks formalism is the need to 
identify which regions of the linear density field are expected 
to collapse to form bound objects. Here one typically gains 
insight from dynamical models that can provide physically- 
motivated thresholds for the linearly extrapolated density 
that is required for collapse to occur. 

One common threshold, based on the spherical collapse 
model (Gunn & Gott, 1972), assumes that haloes can be 
identified in the smoothed linear density field by means of a 
single free parameter: the initial density contrast, Ssc- The 
basic properties of the halo distribution mentioned above 
can then be deduced from the trajectories of the linear den- 
sity field as the spatial scale of the filtering is varied. In this 
case the threshold for collapse at redshift Zc (in an Einstein- 
de Sitter universe) is Ssc ~ 1.686 (1 + Zc)- This result is 
approximately valid for a large range of cosmologies and is 
independent of mass scale and initial density contrast. Eke 
et al. (1996) provide useful analytic formulae for calculating 
the value of 5sc for cosmologies with Q.m + Ha = 1. 

Although the simple spherical model works well in a 
statistical sense its validity is questionable; it is well-known, 
for example, that perturbations in Gaussian random fields 
are inherently triaxial (Doroshkevich, 1970; Bardeen et al., 
1986) . This led to modifications of the spherical model to in 
incorporate the dynamics of ellipsoidal collapse. 

Building upon the classic work by Lynden-Bell (1964) 
and Lin et al. (1965), a number of authors have studied the 
effect of internal shear on the collapse of a homogeneous 
ellipsoid in a uniformly expanding background (e.g., Icke, 
1973; White & Silk, 1979; Peebles, 1980; Watanabe & In- 
agaki, 1991; Lemson, 1993). Eisenstein & Loeb (1995) and 
Bond & Myers (1996) generalized this model to follow the 
collapse of an initially spherical perturbation in the presence 
of external shear. In this case, the sphere is distorted into 
an ellipsoid whose principal axes are parallel to those of the 
external tidal field. One can then ascertain the impact of 
external tidal shear on the properties of collapsed regions. 
For example, Monaco (1997) and Sheth et al. (2001) have 
studied how ellipsoidal collapse modifies the mass function 
of dark matter haloes in the excursion set formalism (see 
also Lee & Shandarin, 1998). 

In the ellipsoidal collapse model of Bond & Myers 
(1996) the collapse time depends explicitly on the local ellip- 
ticity, e, and prolaticity, p, of the tidal field. Less spherical 
regions of the linear tidal field must overcome additional 



tidal stretching and thus require higher overdensities to col- 
lapse. Sheth et al. (2001) showed that the dependence of the 
collapse threshold on e and p could be well-approximated by 
the solution to 
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Here jS and 7 are numerical parameters determined from fits 
to the model results. Sheth et al. (2001) used the results of 
Doroshkevich (1970) to compute the most probable values of 
e and p as a function of smoothing scale. At any ellipticity, 
the most probable prolaticity is Pmp = 0, and for p = they 
find Cmp = o /5/\/Z. Using these results they showed that 
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with /3 = 0.47 and 7 = 0.615. Using their assumptions, eq. 2 
can be expressed as a function of the tidal ellipticity by 
substituting a{M) — ((5/\/5) emp. (Various model-predicted 
collapse boundaries have been compared directly with N- 
body simulations by Robertson et al. (2009).) 

Adding to this, the mass distribution in the vicinity of 
the linear density peaks from which haloes form are inher- 
ently triaxial (Bardeen et al., 1986). This result also applies 
to the Lagrangian regions which collapse onto those peaks 
to form haloes (e.g., Porciani et al., 2002; Ludlow & Por- 
ciani, 2011). This implies that, in the absence of any external 
shear, the internal shear associated with the initially non- 
spherical shape of a collapsing region may impact its col- 
lapse time, and consequently, the average density required 
for complete collapse. These subtleties are often over-looked 
in simple ellipsoidal collapse models, such as those of Bond 
et al. (1991) and Sheth et al. (2001). 

The purpose of the present paper is to investigate how 
simple modifications to the classic ellipsoidal collapse model 
affect the predicted shape of collapse barriers. Our analysis 
focuses on two cosmological simulations of structure forma- 
tion, the results of which we compare directly with predic- 
tions from the model. The rest of the paper is organized 
as follows. In Section 2 we introduce our simulations and 
briefiy describe our main analysis techniques. Section 3 de- 
scribes some of the basic characteristics of proto-haloes in 
our simulations, focusing on their inertial and tidal tensors. 
We use these results in Section 4 to motivate a closer look 
at the ellipsoidal collapse model, and in particular, how the 
predicted collapse barriers depend on the assumed shape of 
the collapsing region. The predictions are compared directly 
with the results of two cosmological simulations in Section 5. 
Finally, we end with a brief discussion in Section 6 before 
summarizing our main results in Section 7. 



2 NUMERICAL ISSUES 

We present here, for completeness, a brief summary of our 
numerical simulations and analysis techniques. More details 
can be found in Pillepich et al. (2010) and Ludlow & Porciani 
(2011). 
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2.1 The simulations 

We focus our analysis on two high-resolution cosmological 
simulations of structure formation in the standard LCDM 
cosmology. Our runs adopted cosmological parameters con- 
sistent with the WMAP five-year data release (Komatsu 
et al., 2009). These parameters are as follows: Q.m = 0.279, 
^A = 1 - ^M = 0.721, as = 0.817, Us = 0.96, and a Hubble 
constant Hq = H{z = 0) = 73 km s^^ Mpc. 

Both runs followed the evolution of the dark matter 
component in comoving coordinates using a lean version 
of the simulation code Gadget (Springel, 2005). The lin- 
ear density field was sampled using 1024^ equal mass par- 
ticles in a periodic box with side-lengths equal to l^yox = 
150 h~^ Mpc and 1200 h~^ Mpc. Initial conditions were 
generated using the the Zel'dovich approximation, adopting 
a starting redshift Zs = 70 for the 150 h^^ Mpc box run, and 
Za = 50 for the 1200 h~^ Mpc box (a thorough discussion of 
our initial conditions can be found in Pillepich et al. (2010)). 
For these choices the particles masses are 2.433xl0*/i^^ 
Mq and 1.246x10" h'^ M© in the 150 ft" ^ Mpc and 
1200 h~^ Mpc boxes respectively. The corresponding soft- 
ening lengths are 3 h^^ kpc and 20 h^^ kpc, which are kept 
fixed in comoving coordinates during the simulation. 



2.2 The analysis 

We study some of the basic properties of the Lagrangian 
regions that eventually collapse to form haloes. We refer to 
these as proto-haloes for convenience. Dark matter haloes 
were identified at z = using a friends-of-friends (FOF) 
halo finder with linking-length I — 0.2 times the mean inter- 
particle separation; the "proto-haloes" are identified in the 
simulation initial conditions by tracing backwards in time 
the subset of particles that belong to these FOF haloes. 
The proto-halo mass is defined as the total mass associated 
with a given FOF halo. We have verified that our results are 
insensitive to our adopted halo definition by repeating the 
analysis with dark matter halos identified with a spherical- 
overdensity halo finder. These halos were identified at a den- 
sity contrast level of 200 x (p) . 

We quantify proto-halo shapes using the following de- 
scription of their mass distribution (Cole & Lacey, 1996; 
Bailin & Steinmetz, 2005; Hopkins et al., 2005; Shaw et al., 
2006; Rett et al., 2007): 
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where rk,i is the distance from the center of mass to parti- 
cle i, and the factor of l/r| is present to prevent particles 
at large distances from dominating the shape estimates. The 
center of mass is defined, at the initial redshift of the simula- 
tion, as the geometric center of all FOF particles associated 
with the z = halo. 

Eq. 3 can be diagonalized and the axis lengths, ai > 
0-2 > ffl3, determined from the square roots of the eigenval- 
ues; the eigenvectors coincide with the principal axis frame 
of the system. Once the axis lengths Ui have been found the 
shapes of the proto-haloes can be characterized in terms of 
their ratios: a-ijax measures the intermediate-to-major axis 
ratio, and as,/a\ the minor-to-major. 

In addition to shapes, we also measure properties of 



the linear tidal field at the geometric center of each proto- 
halo. The linear tidal field is calculated from the deformation 
tensor: 

(4) 



'^ dxidxj 

Here "l?(x) is the peculiar gravitational potential at position 
X, which is related to the density contrast through Poisson's 
equation: V^<I> = 5. The ordered eigenvalues of O, Ai > A2 > 
A3, determine the geometry of the linear density field, and 
their signs indicate whether the fiow around a given point 
is of compression (-I-) or dilation (— ). 

The eigenvalues of T> can also be expressed in terms of 
the tidal ellipticity, e, and prolateness, p, defined 

A1-A3 
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and 
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where 5 — Ai -|- A2 + A3. For 5 > 0, we have e > and 
— e < p < e. 

In what follows we will measure the eigenvalues of the 
linear deformation tensor at proto-halo barycenters after av- 
eraging the relevant quantities over a spherical Lagrangian 
volume enclosing the halo mass. 



3 PROTO-HALOES IN COSMOLOGICAL 
SIMULATIONS 

In Fig 1 we plot the mass dependence of the axis ratios, 
0,2/0.1 and 03/ a\, for all proto-haloes in both of our simu- 
lations. Colored points show the median trends; boxes and 
whiskers highlight the lO"', 25*'', 75 "^ and 90"^ percentiles 
of the distribution. Results for haloes in our 150 h~^ Mpc 
box are plotted in blue, and red points are used those in our 
1200 fe"^ Mpc box. 

There is a clear mass dependence to the mean proto- 
halo shape. Judging their sphericity by the value of 03/01, 
we find that, on average, (03/01) varies from ~ 0.4 for haloes 
with Mpop ~ lO^^/i~^M0 to ~ 0.65 for rich galaxy clusters, 
albeit with considerable scatter an any given mass scale. It 
is important to note that spherical proto-haloes are quite 
rare in our sample. In our 150 h~^ Mpc box only ~ 0.4% 
of proto-haloes containing more than 500 particles have 
03/01 > 0.8. For our larger box, which probes higher mass 
scales, ~ 2.9% of (Npop > 500) haloes have 03/01 > 0.8. 
Thus, the Lagrangian regions subject to the collapse con- 
straint are highly non-spherical. 

For haloes containing more than 500 particles (indicated 
by vertical arrows in Figure 1) the median trends for each 
run agree quite well and can be approximated by: 



— = Q X logm(MFOF/[lO'°/i-'M0]) + So, 

Ol 



(7) 



as shown by the dotted lines in Fig 1*. The best-fit pa- 
rameters, Q, and So are listed in table 3. Note also that 



Note that the best-fit parameters in cq 7 vary only slightly 
when halos arc identified as spherical-overdensitics with density 
contrast 200 X (p). Best-fit parameters for both FOF and SO halos 
are given in Table 3. 
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axis ratio Halo Dcf. Q 



So 



a^j Ci\ 


FOF 


0.0286 


0.666 


a^/ai 


FOF 


0.0618 


0.345 


02/01 


SO 


0.0277 


0.671 


as/ai 


SO 


0.0629 


0.336 



Table 1. Best-fit parameters describing the average axis ratios, 
02/01 and 03/01, for dark matter proto-haloes as a function of 
mass. Results for both FOF and SO halos are provided. 



proto-haloes are nearly triaxial, as indicated by the thick 
black lines in the lower panel of Fig 1. This line shows the 
ratio as/ai expected for "perfectly triaxial" ellipsoids (i.e. 
those with a2/ai = 0,3/0.2) given the median values of 02/01 
in the upper panel. In what follows we shall only consider 
haloes with A^fof > 500, unless explicitly stated otherwise. 

What determines the shape of a proto-halo? Sev- 
eral studies have assumed that proto-halo boundaries can 
be identified with a fixed-level iso-density contour drawn 
around local density maxima in the linear density field (e.g., 
Bardeen et al., 1986; Hoffman, 1986, 1988; Heavens & Pea- 
cock, 1988; Catelan & Theuns, 1996). However, as discussed 
by Porciani et al. (2002), this scenario is not true for proto- 
haloes, in general; the boundaries of proto-haloes are not 
fully determined by the self-gravity of the local mass distri- 
bution, but rather by the tidal field generated by the external 
mass distribution. 

In Figure 2 we show an example of this. Here we plot 
the linearly extrapolated density field, 1+5 centered on a 
2.99 X lO^^/i"^M0 proto-halo (red dots) identified in the 
initial conditions of our 150 h~^ Mpc box. Isodensity con- 
tours corresponding to fixed density contrasts of (5 = 1.0 and 
5 — 1.686 are also shown. The iso-density contours enclosing 
the peak around which the halo collapses provides an intu- 
itive sense of the (projected) spatial distribution of dark 
matter surrounding the peak: the peak is extended along 
a high density ridge connecting two massive structures and 
separating two low density voids. The long axis of the proto- 
halo, however, lies transverse to that of the density ridge in 
which it forms. It is easy to understand why. Material flow- 
ing out of the voids leads to maximal compression along 
an axis roughly perpendicular to the density ridge, whereas 
tides from the surrounding large scale structure cause di- 
lation along the ridge. The net result is a push-and-puU of 
material that sets the shape of the proto-halo such that its 
long axis coincides with the direction of maximum compres- 
sion, or infall, and its shortest axis to direction of maximum 
dilation. The grey arrows in Figure 2 show the linear velocity 
field in the vicinity of the proto-halo, and help to clarify the 
impact of large-scale tides on gravitational collapse around 
density peaks. 

Similar results were discussed previously by Porciani 
et al. (2002), who found a very strong alignment between the 
eigenvectors of proto-halo's inertia tensors and those of the 
tidal shear at the proto-halo's center. Their results, based on 
one cosmological N-body simulation, indicate that the the 
geometry of the linear tidal field has a substantial impact 
on the geometry of Lagrangian regions which are destined 
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Figure 1. Mass dependence of the axis ratios, a2/a\ and 03 /oi, 
computed for proto-haloes in our simulation initial conditions. 
Points show the median trends; boxes and whiskers the lO***, 25*'', 
75*'', and gO*** percentiles of the distribution in each mass bin. 
Blue symbols correspond to our 150 h~^ Mpc box simulation; red 
to our 1200 /i~^ Mpc box. Vertical arrows indicate the lower 500- 
particle limit imposed in our halo sample. The dotted lines show 
the best-fit (eq. 7) to the combined median trends for mass bins 
corresponding to haloes with more than 500 particles. The thick 
black lines in the lower panel show the minor-to-major axis ratios 
expected for "perfectly triaxial ellipsoids": 02 /oi = 03/02. 
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Figure 2. Linear density and velocity fields in the vicinity of a proto-halo identified in our 150 h~^ Mpc box simulation. The proto-halo 
is shown using red dots; thick contours highlight fixed thresholds, 5 = 1.0 and S = 1.686, of the linear density field, which has been 
smoothed on the halo mass scale and linearly extrapolated to 2 = 0. The halo forms in the vicinity of a density peak of the same 
characteristic mass, yet the mass distribution of the proto-halo and it associated density peak are very different; the geometry of the 
peak is determined only by the density threshold, whereas the proto-halo shape and orientation are determined by the external tidal 
field generated by the surrounding large scale structure. 



to form haloes. Whether or not the inertia tensors of den- 
sity peaks (above a given density contrast) correlate with 
their local tidal field remains a pending issue which holds 
particular relevance for theories of galaxy formation. 

From a dynamical perspective the in and outflow of ma- 
terial around proto-halo centers can be characterized by the 
signs of the eigenvalues of the deformation tensor measured 
at their centers of mass. These indicate whether the initial 
flows along the principle axes are of compression {+) or di- 
lation (— ). 

If proto-halo boundaries are set by the tidal field then 
both the relative magnitudes and the signs of the A^.s de- 
termine their shape. In Figure 3 we plot the ratios ai/uj 
versus \i/\j for proto-haloes in our simulations, showing 
separately two different mass bins. The AiS are measured at 
each proto-halo's geometric center after smoothing the lin- 
ear tidal field on the halo mass scale with a spherical tophat 
filter. In our 150 h~^ Mpc box we select haloes in the range 
10^^ < MFOF/(/i~^Af0) < 5x10^^; the contours show the 50 
and 75 percentiles of their distribution. In our 1200 h~^ Mpc 



box 2686 haloes are selected from a higher mass bin span- 
ning 5 X lO"/i"^M0 to 2 X lO^^h-^MQ. The medians of 
the joint distributions for these mass bins are shown as blue 
circles and red triangles, respectively. The Spearman rank 
correlation coefficients for each of these samples are pro- 
vided in the panels. In the high mass bin all proto-haloes 
have Ai > A2 > A3 > 0; in the lower mass bin Ai > A2 > 
for all haloes, but A3 > for only 41. 9%"^. 

The trends in Figure 3 can be interpreted in terms of 
the V — I correlation discussed at length by Porciani et al. 
(2002) and Lee & Pen (2000). Consider each of the columns, 
in turn. In the left-hand panels we plot the various axis ratios 
ai/uj versus A2/A1. For A2 ~ Ai roughly equal compression 



T In general, we find that essentially all proto-haloes with 
NpOF > 100 have Ai > A2 > 0, but that there is a strong mass 
dependence to the fraction with A3 > 0. Below 10^^/i~^Mq, for 
example, ^ 87% of proto-haloes have A3 < 0; above lO^*/i~^M0 
^ 97% have A3 > 0; A3 > for 50% of haloes at MpoF ~ M*, 
the characteristic mass for collapse at 2 = 0. 
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occurs along the major and intermediate axes of the tidal 
field. In this case we might expect the proto-halo to be ap- 
proximately spherical (in the case that A3 ^ A2), or oblate 
(if A3 ^ 0). The Massive haloes (red curves) with A3 > are, 
on average, more spherical than the lower mass ones: the av- 
erage value of (02/01) is ~ 0.80, and for (03/01) ~ 0.65, with 
very little mass-dependence to the mean. In the lower mass 
bin, however, haloes have —0.5 ^ A2/A1 ;$, 0.5 which leads 
to weak trends between the axis ratios and the ratio A2/A1. 
These trends can be interpreted as follows. For A2 ~ Ai 
haloes are slightly oblate; 03/02 decreases as A2 — > Ai, while 
02/ai increases. As A2 — >■ 0, however, the systems then be- 
come more prolate, since Ai > A2 ~ 0, which guarantees 
A3 ^ 0. In this case, 02/01 decreases while 03/02 increases 
(more prolate configurations). 

Next we consider the far right panels, which plot the 
ratios Oi/oj against A3/A2. Again certain correlations are 
seen amongst these parameters. For A3 ~ A2 proto-haloes 
are close to spherical (if Ai ~ A2 ~ A3) or prolate (if Ai > 
A2 ~ A3). For prolate systems 03/02 increases as A3 — )■ A2; 
for spherical systems both 02/01 and 03/02 increase. For 
small or negative values of A3 we see a sharp decrease in the 
ratios 03/01 and 03 /02 as the systems become less spherical, 
and less prolate, respectively. 

Finally we consider the middle panels, which plot o^/oj 
verses A3/A1; these show the strongest and least mass- 
dependent correlations of all. The trends can be under- 
stood as follows. In cases where A3/A1 is large and posi- 
tive the compression exerted by external tidal forces is most 
isotropic and c — >■ 6 — >■ o (recall that in the spherical limit 
Ai,2,3 = S/3). As A3 becomes large and negative, however, 
there is substantial dilation along the minor axis of the tidal 
field. Since A2 > for all of our haloes we expect to see a 
sharp decrease in the axis ratios 03 /oi and 03 /02 due to the 
tidal dilation along the minor axis of the proto-halo. 

It should be emphasized that these correlations are sub- 
ject to the fact that the proto-haloes are destined to collapse 
to form a single object hy z = 0. The Lagrangian geometries 
of other present-day structures, such as filament, sheets or 
voids may differ substantially from these. 



4 NON-LINEAR DYNAMICS: ELLIPSOIDAL 
COLLAPSE OF A TRIAXIAL 
OVERDENSITY 

Based on these results, we consider here modifications to the 
classic ellipsoidal collapse model of Bond & Myers (1996) in 
order to follow the collapse of an initially triaxial (rather 
than spherical) perturbation in the presence of an external 
tidal field. In principle, this situation is much more com- 
plicated since the external shear need not be aligned with 
the principal axes of the initial ellipsoid. However, in their 
systematic study of proto-halo shapes Porciani et al. (2002) 
found a nearly perfect alignment between the principal axes 
of the total shear and the inertia tensor of proto-haloes in 
numerical simulations (see also Lee & Pen, 2001; Lee et al., 
2009). This suggests that either the internal shear domi- 
nates over the external, or that the shapes of proto-haloes 
are dictated by the external shear. 

In our model we consider a uniform density pertur- 
bation on top of a fiat Friedmann-Robertson- Walker back- 
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Figure 3. Axis ratios of proto-haloes plotted versus the ra- 
tios of the eigenvalues of the deformation tensor. Blue curves 
(with circles) show the medians for ^ 10* haloes selected from 
our 150 h~^ Mpc box which have masses in the range 10^^ < 
MpQp/{h~^MQ) < 5 X lO'^^. Red curves correspond to haloes 
with 5 X IQl" < Mpof/(/j"^M0) < 2 x lO^^s; there are 2686 
of these. The Spearman rank correlation coefficients for the joint 
distributions are given in each panel. The contours enclose the 50 
and 75 percentiles of the distributions for haloes in the low mass 
bin. In all cases, the A^s are measured at the proto-halo barycen- 
ter after smoothing the tidal field with a tophat filter containing 
the halo mass. 



ground whose energy content is dominated by the matter 
density, pm, and a cosmological constant, A. We model the 
perturbation as an homogeneous ellipsoid with semi-axes of 
physical length ai (i = 1,2,3) and density contrast 5, and 
adopt Cartesian coordinates, Xi, measured along the prin- 
cipal axes of the ellipsoid. In this coordinate system, the 
internal gravitational potential of the ellipsoid can be writ- 
ten 



$oll(x) = TtG (1 + 5) pni ^ hxi 



(8) 



where the coefficients bi are given by Carlson's elliptic inte- 
grals (Kellogg, 1929; Chandrasekhar, 1969): 



oi 02 03 



dr 



(9) 



° {r + aj)\{{i 



2x1/2 



J = l 



For a sphere 61 = &2 = &3 = 2/3 and, in general, ^ &i = 2 
as a consequence of Poisson's equation. The potential gen- 
erated by the uniform distribution of matter outside the 
perturbation can be written 



<I>(x) = -TvGpm y^ 



Ac^' 



(10) 
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Figure 4. Evolution of the axis lengths for collapsing ellipsoids. 
Left panels assume the initial overdensity is a uniform sphere; 
right-hand panels adopt an initially triaxial overdensity with axis 
ratios a-^jax = 03/02 = 0.8. The evolution of the axis lengths 
is determined by the total tidal field, (e,p, (5), as shown in each 
panel. For the initially triaxial perturbation we modify the exter- 
nal tidal field to account for the initially non-zero internal com- 
ponent, such that the total tidal field is the same in the left and 
right-hand panels. 



Here fii — bi — 2/3 and c is the speed of light. 

In a reahstic cosmological setting the perturbation will 
be immersed in a randomly fluctuating density field and will 
be subject to external tidal forces. The net effect of these 
tides can be approximated by expanding the external poten- 
tial in multipoles and focusing on the quadrupole term. (The 
dipole term generates a translation of the ellipsoid's cen- 
ter of mass but does not alter its shape; contributions from 
higher-order moments are generally subdominant (Quinn & 
Binney, 1992)). In a coordinate system which is parallel to 



the principal axes of the quadrupole, the tidal gravitational 
potential is then $tid(y) = J2i '^f"*'^?. with X), A,'""'' = 0. 
Assuming that the tidal shear is parallel to the princi- 
pal axes of the density perturbation (which naturally arises 
when one starts with a sheared spherical perturbation), we 
can write the equation of motion for the semi-axis lengths 
of the ellipsoid as 



at — 






(cxt) 



(11) 



where dots denote time derivatives. 

Since the potential is quadratic, and the acceleration is 
linear in the coordinates, all ellipsoids with the same axes 
aspect ratios evolve self-similarly, and the density pertur- 
bation remains homogeneous. Given the comoving coordi- 
nates of the axes extrema, qi, the evolution thus preserves 
fli a2 as (1 + 5) = qi q2 qs a^ = constant, and so 

5^Sl^l3La^_l, (12) 

ai a2 a-i 

Initial conditions can be set at some early time, to, using 
the Zel'dovich approximation: 

ai{to) = ?, [l-A>(fo)] a{to) , (13) 

di{to) = H{to)ai(to)-qiHD{to)Xr{to)a{to) . (14) 

Here H = a/ a is the Hubble parameter; D the linear growth 
factor; Ho = D/D, and \i(t) oc D{t) are the eigenvalues of 
the linear deformation tensor. In order to match the last 
expressions with eq. 11, we also derive the eigenvalues of 
the deformation tensor at io from the total gravitational 
potential: 



Ar'(io) = '5(to) 



1 Pijta] 
3 2 



+ A 



(cxt) 



(io) 



(15) 



Only the components of this expression that scale as D{t) 
will contribute to the Ai terms appearing in the Zel'dovich 
approximation. 

The early evolution of the internal gravitational poten- 
tial of the ellipsoid can be determined by combining eq.s 9 
and 13 and then performing a Taylor expansion in powers 
of D(t): 

A(i) = Pf^ + A/3, D{t) + 0[D\t)] . (16) 

Here j3i is evaluated using eq. 9 with ai — qi. The second 
term in the expansion is given by 



3 
A/3i D{t) = ^ Uij \j 



j=i 



where 



ctjj = gi 172 q-i 



and, for j 7^ k, 



Q-jk = —q\ q2 qs 



(2g: 



t) dr 



(17) 



(18) 



ilj+rrl[{q'+r: 



1/2 



r dr 



{<l^+r){ql+r)l[{q? + r: 



(19) 



1/2 



When all the gi's are equal we have Oij — 8/15 when i = j 
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Figure 5. Comparison of tiie collapse boundaries predicted by tlie ellipsoidal collapse model for various initial perturbation shapes. 
The different panels show separately the linearly extrapolated over-density required for collapse along the primary axis (top-left), the 
intermediate axis (top-right), and the minor axis (bottom-left) of the initial perturbation. The bottom-right panel shows the barrier 
height required for all three axes to collapse. The shaded and hatched regions show the effect of varying the prolaticity of the tidal field 
over the range — e/2 < p < e/2. Orange, blue and white regions show results for different initial perturbation shapes: a2/ai = 03/02 = 1, 
0,2/0.1 = 13/02 = 0.8, and 02/01 = 03/02 = 0.6, respectively. In all cases, the major axis of the tidal field is assumed to be aligned with 
the major axis of the initial perturbation; the corresponding intermediate and minor axes are also aligned. 



and —4/15 when i 7^ j. Note that this reduces to the re- 
sult of Bond & Myers (1996), A ~ Al3i D(t) = (4/5) ti for 
qi — q2 ~ qs- Thus, the term P5/2 in eq. 15 is of second 
order in D{t) and does not contribute to the linear shear 
tensor, which is fully generated by the external tides. How- 
ever, contrary to the case discussed by Bond & Myers (1996), 
for initially triaxial perturbations the linear velocity shear 
is composed of two terms 



.4- 



PI 



(0) 



-+Ar' 



(20) 



For the typical axes ratios of proto-haloes in our simulations 
{q2/qi ^ q3/q2 ^ 0.8) we find 



Pi"' ~ (-0.170,-0.011,0.182), 



(21) 



showing that internal shear slows down collapse along the 
long axis and speeds it up along the short axis. 

Solving eq. 11 requires knowledge of the time evolution 
of the external tidal field, Af'"'. Bond & Myers (1996) pro- 
posed two limiting approximations. In one case, the external 
structure is assumed to grow independently of the ellipsoid 
and the shear scales as the linear growth factor. 



A 



(cxt) 



it) 



Djt) 
D{to) 



(cxt) 



(io) 



(22) 



This approximation is accurate when the shear is dominated 
by the contribution from large scales. Alternatively, one can 
assume that the external shear is dominated by the contribu- 
tion from the regions immediately surrounding the ellipsoid, 
and is completely driven by the shape of the ellipsoid itself. 
In this case, Xf^ (t) — (5/4) /3i (note that this reduces to 
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eq. 22 at early times). The relation above only holds for 
an initially spherical perturbation but can be generalized to 
triaxial ones using 



Ar'w=E"^'(/3^-^r) 



(23) 



where the a~ are the elements of the inverse matrix of aij . 
A hybrid model which interpolates between these asymp- 
totic regimes has been recently proposed by Angrick & 
Bartelmann (2010). In this case one follows the non- linear 
growth proposed by Bond & Myers (1996) until the corre- 
sponding axis turns around while the linear growth is used 
afterwards. For simplicity, we consider only the linear evo- 
lution of the external tides in this work. 

The ellipsoidal model for halo collapse is not fully self 
consistent. For example, the presence of the ellipsoidal per- 
turbation will produce inhomogeneities in the background 
(the potential is not quadratic outside the perturbation) 
which will, in turn, perturb the density distribution within 
it. This effect is usually neglected by assuming that both 
the perturbation and the background remain homogeneous 
at all times (e.g., Icke, 1973; White & Silk, 1979; Eisen- 
stein & Loeb, 1995). Furthermore, external tides cannot be 
generated by a smooth background, which should be inter- 
preted as an effective mean field. The presence of substruc- 
ture (which are not considered in the ellipsoidal model) can 
also modify (and in particular slow down) the collapse of 
a perturbation (e.g., Davis & Peebles, 1977). Finally, the 
equations of motion do not account for the physics that re- 
sults in the virialization of the collapsing structure. These 
relaxation processes are crudely mimicked by "freezing" the 
collapse along each axis once a critical radius aoq,i — f aqi 
is reached. Bond & Myers (1996) suggested using the ra- 
dial freeze-out factor / = 0.177, which reproduces the virial 
density contrast of ~ 178 characterizing the collapse of a 
spherical perturbation in an Einsten-de Sitter universe. (In- 
troducing this radial freeze-out factor also modifies slightly 
the characteristic density required for collapse; in the spher- 
ical limit Sc ~ 1.686 for collapse to an infinite density point, 
<5c ^ 1-609 for / = 0.177.) More recently, Angrick & Bartel- 
mann (2010) have used the tensor virial theorem to predict 
more realistic values of /. However, the last phase of col- 
lapse is generally rapid, and altering the value of / does not 
strongly affect the estimated epoch of virialization. For ease 
of comparison with previous work, all applications of this 
model presented below have assumed an Einstein-de Sitter 
background cosmology, and a freeze-out factor / — 0.177. 

Figure 4 plots the trajectories of the three axes lengths, 
in physical units, for the ellipsoidal collapse model of Bond 
& Myers (1996) (left panels) and for the modified model 
described above (right). We start our numerical calculations 
at to = 0.250 Myr, at which point 5{to) = 1.686 sets the 
overdensity based on the spherical collapse criterion. The 
top two panels show the evolution of an initial spheroidal 
(left) and ellipsoidal perturbation with axis lengths 02/01 — 
0.3/0,2 = 0.8 (right) when no external tides are present. The 
lower panels show, for the same initial set-up, the effect of 
varying the external tidal field on the evolution of the axis 
lengths. In all cases the external tides are imposed at Zi so 
that the total tidal field acting on the overdensity is the same 
in both models. 



Figure 4 highlights a number of things worth noting. 
First, for a given e, the Bond & Myers (1996) model pre- 
dicts the formation of pancake-like structures when p > 0; 
for p < filamentary structures form. Thus, for a given 
(e,p) the ellipsoidal collapse model predicts the axis ratios 
of the object at all times, and in particular, at the final 
time. These results have been used recently by Rossi et al. 
(2011) to model the shape distribution of collapsed struc- 
tures within the excursion set theory. Notice, however, that 
relaxing the assumption of spherical symmetry at Zi sub- 
stantially changes the evolution of the axis lengths, even 
when the total tidal field is chosen to be identical. In par- 
ticular, the ordering of the axes lengths can change; what 
is initially the primary axis, for example, can swap between 
intermediate and minor, and then back to primary again. In 
addition, the turnaround and freeze-out times also differ in 
the initially triaxial model, a result which, as we show below, 
has important implications for modeling collapse boundaries 
for halo formation. 

In Figure 5 we illustrate the role played by the initial 
shape of the ellipsoidal perturbation in setting the collapse 
thresholds. Orange, blue and white regions show the results 
for 02/01 = 03/02 = 1, 0.8 and 0.6, respectively; the shad- 
ing shows the effect of varying the tidal prolaticity over the 
range — e/2 < p < e/2. Each panel plots separately the 
linearly extrapolated density contrast required for collapse 
along each of the three principle axes of the initial perturba- 
tion (the first axis corresponds to the major axis of inertia; 
the third to the minor axis). In the lower-right panel we 
plot the "collapse boundary", defined as the linear density 
contrast required for all three axes of the initial ellipsoid to 
collapse. 

In the absence of external tidal forces an initially non- 
spherical perturbation requires higher density contrasts for 
collapse than the spherical case. This is because the internal 
shear acting along the long axis of the ellipsoid always slows 
collapse, and external tides are required in order to balance 
the internal forces . Because p = corresponds to A2 = 0, 
collapse along the intermediate axis is not strongly affected 
by tides; only the (sub-dominant) internal shear induced by 
the aspherical nature of the initial perturbation matters for 
this. Along the minor axis of the initial ellipsoid collapse 
barriers are always lower than the spherical case. Tthis is 
because the internal shear acting along this axis always aids 
in collapse, and very large values of tidal ellipticity (cor- 
responding to very negative values of A3) are required to 
suppress this. 

The collapse barriers inferred from the ellipsoidal col- 
lapse model depend strongly on the assumed initial shape of 
the collapsing ellipsoid. In the presence of the same tidal field 
the time required for a given axis to collapse is sensitive to 
its initial displacement from the perturbation center. Conse- 
quently, the linearly extrapolated over-density required for 
collapse can vary significantly for different initial shapes. 



■f Recall that the primary axis of the perturbation is assumed 
to align perfectly with the major axis of the tidal field, so ' 
ellipticities correspond to greater compression along this axis. 
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Figure 6. Smoothed linear overdensities, extrapolated to 2: = 0, 
plotted as a function of the ellipticity of the tidal shear for proto- 
haloes containing more than 500 particles. Both quantities are 
measured at the proto-halo center of mass and are averaged over 
a spherical Lagrangian region containing a mass MpoF- Colored 
circles show the predicted collapse thresholds derived from the 
ellipsoidal collapse model for several different initial perturbation 
shapes. In each case, the initial perturbation is assumed to be 
"perfectly triaxial" (i.e. 02/01 = 03/02) and we have set p = 
for simplicity. The solid lines (of the same color) show the best-fit 
collapse boundary given by eq. 24; best-fit parameters are given 
in Table 2; residuals from the best-fits are shown in the bottom 
panel. For comparison, we also show the Sheth et al. (2001) barrier 
using a dashed line. 



5 COLLAPSE THRESHOLDS FROM 
ELLIPSOIDAL COLLAPSE 

In this Section we wish to give a more quantitative assess- 
ment of the effects previously discussed. In what follows, 
we will consider the collapse barriers inferred from the el- 
lipsoidal collapse model and express them as a function of 
the ellipticity of the tidal field, e. We will compare these 
directly with the actual overdensities and tidal ellipticities 
measured at proto-halo centers of mass, with both quantities 
averaged over the appropriate Lagrangian volume. Unless 
explicit stated otherwise, we shall follow Sheth et al. (2001) 
and assume p — when calculating barrier heights from the 
ellipsoidal collapse model. 



Figure 7. Same as Figure 6 but for haloes in four narrow mass 
bins. Contours enclose the 50 and 75 percentiles of the measured 
collapse boundaries for haloes in each mass bin. In each panel the 
dot-dashed curve shows the ellipsoidal collapse boundary for an 
initially spherical perturbation acted on by external tidal forces 
with p = and ellipticity e. The solid curves show the barrier 
heights for an initially triaxial perturbation with axis ratios con- 
sistent with the average axis ratio of haloes in each mass bin. 
The total tidal field (internal plus external) is the same in both 
cases. For comparison, we also show the fitting formula of Sheth 
et al. (2001) as a dashed line, converted to the most probable halo 
ellipticities via Cmp = (j/S/\/5. 
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Table 2. Best-fit parameters. A, Bi, and B2, for the ellipsoidal 
collapse barriers (eq. 24) shown in Figure 6. The model results 
are based on the collapse of "perfectly triaxial" ellipsoids (i.e. 
02/01 = 03/02) embedded in a tidal field for which p = 0. 



5.1 Collapse barriers from the ellipsoidal collapse 
model 

Our adaptation of the ellipsoidal collapse model of Bond & 
Myers (1996) yields a barrier height that is both a function 
of the statistics of the tidal field at any given point, (e,p), 
and also of the sphericity of the collapsing region. To reduce 
the number of degrees of freedom in this model we make the 
following simplifying assumptions: a) CDM proto-haloes can 



be modelled as "perfectly triaxial" ellipsoids (i.e. 02/01 = 
03/01), and b) the prolaticity of the tidal field takes on the 
mean value at any mass scale, i.e. p = 0. 

Figure 6 plots the various collapse boundaries, S{e), in- 
ferred from our ellipsoidal collapse model for several initial 
perturbation shapes. Different colored points correspond to 
initial overdensities with different axis ratios, which are in- 
dicated on the plot. The data are well reproduced by the 
simple function 
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A [exp(-Bi e) +exp(B2 e)]^ 



(24) 



as can be seen by the solid curves in Figure 6, whose colors 
are chosen to match the model predictions. These curves 
differ markedly from the Sheth et al. (2001) barrier, which 
is shown with a dashed black line. The best-fit parameters 
vary systematically with the initial shape of the ellipsoidal, 
and are listed for each of the plotted models in Table 2. 

In the lower panel of Figure 6 we plot the residuals 
from the best-fit to eq. 24. The residuals are typically small, 
^ 5% over a wide range of ellipticities. Note, however, that 
the model predictions have a much sharper upturn than the 
best-fit profile; beyond some characteristic ellipticity the fits 
diverge from the model predictions. However, as we have al- 
ready seen, the average minor-to- major axis ratio varies from 
(as/ai) ~ 0.4 at halo masses ~ lO^^/i^^M0 to ~ 0.65 on the 
scale of rich galaxy clusters. The mass dependence of the av- 
erage value of (02/01) is weaker, ranging from ~ 0.7 to ~ 0.8 
over the same mass range. On average, the most appropriate 
barriers to compare to the simulation results are, therefore, 
the yellow and green curves in Figure 6. For these, the bar- 
rier heights predicted by the ellipsoidal model are accurate 
to within ~ 5% for ellipticities between roughly 0.65 and 0.8; 
in our 150 h~^ Mpc-box simulation only ~ 6.8% of haloes 
have e > 0.65; 2.7% have e > 0.8. The typical mass of these 
objects is low, with an average of ~ 32 x lO^^/i'^M© . For the 
higher mass scales probed by our 1200 h^^ Mpc run, there 
are no haloes with e > 0.65. Our fits are therefore accurate 
over the range of ellipticities and halo masses probed by our 
simulations. We discuss this issue further in the following 
subsections. 

The shaded area and contours show the distribution of 
{5h, eh) for all proto-haloes with A''fof > 500 in both of our 
simulations. Clearly the very large scatter in the average 
linear overdensity at dark matter proto-halo centers is not 
captured by the standard version of the ellipsoidal collapse 
model. Considering the large variation in barrier heights con- 
sistent with the ellipsoidal collapse model it is interesting to 
investigate the extent to which the variation is induced by 
initial shapes of proto-haloes. We turn to this next. 



5.2 Mass dependence of barrier heights 

Figure 7 plots the dependence of the Barrier height on 
tidal ellipticity for haloes selected to lie in several nar- 
row mass bins spanning the range Mfof ~ IQ^^h^^VlQ 
to ~ 6 X lQ^*'h^^yiQ. Red contours enclose the 50 and 75 
percentiles of the distributions for our simulated halo data; 
the curves show various theoretical estimates of the bar- 
rier height. The dashed curves show the fit given by Sheth 
et al. (2001) assuming e — emp ~ (J /5/^\ the dot-dashed 
curves shows the results of the "standard" ellipsoidal col- 
lapse model (with p = and qi = q2 = Hi)', the solid curves 
shows the barrier height computed using the average (92/^1) 
and (qa/qi) for haloes in each bin. The percentage of halos 
whose measured overdensities exceed these barriers are pro- 
vided next to each curve. 

Clearly the minimum barrier height for haloes of any 
given mass is more accurately described by the ellipsoidal 
collapse model once the initial perturbation shape is cho- 
sen to match the mean shape in each mass bin. At the 
lowest masses, where proto-haloes can have the most flat- 



tened geometries, only ~ 49.6% of proto-haloes lie above 
the standard ellipsoidal collapse barrier, whereas ~93.8% 
lie above the collapse barrier inferred for proto-haloes with 
02/01 ~ 0.70 and 03/ai = 0.45, the mean values for this 
mass. Towards higher masses, where proto-haloes are intrin- 
sically more spherical, we find ~ 73.2% lie above the dot- 
dashed curve (56.0% above the Sheth et al. (2001) barrier) 
and ~ 95.5% above the barrier when modelled to account 
for the average initial shape of the Lagrangian proto-halos. 
Note, however, the large scatter in (e^, Sh) at any given 
mass scale. Within the standard ellipsoidal collapse model 
described by Bond et al. (1991), the scatter in barrier heights 
at any mass and ellipticity is completely determined by the 
prolaticity of the tidal field in the vicinity of the proto-halo. 
The observed scatter in Sh at fixed e^ shown in Figures 6 
and 7 can not be accounted for by variation in prolaticity 
alone. By comparing the scatter amongst the red points in 
Figure 7 to the grey shaded regions in the lower-right panel 
of Figure 5, it is clear that changes in prolaticity alone can- 
not account for the variety of measured barrier heights at 
fixed halo mass. In the next section we explore whether the 
variation in proto-halo shape can account for some of this 
scatter. 



5.3 Shape dependence of barrier heights 

Some insight into the origin of scatter in barrier heights 
can be gained by considering how the measured proto-halo 
boundaries vary as a function of proto-halo shape at fixed 
mass. In Figure 8 we plot the joint distribution of {5h,eh) 
for proto-haloes in our 150 h^^ Mpc box simulation that 
lie in the mass range 24.3 < MFOF/(lO^°h"^M0) < 48.7 
(the same proto-halo sample plotted in the upper-left panel 
of Figure 7). These proto-haloes have been split into four 
separate subsamples based on the value of their axis ratios. 
In order to make a more meaningful comparison with the 
predicted collapse barrier, which we have computed for a 
"perfectly triaxial" ellipsoid, we also select proto-haloes to 
have 02/ai « 03/02 spanning a narrow range. Proto-haloes 
in the upper-left panel are the most spherical in this mass 
bin, and 03/01 > 0.8; in the upper-right panel all proto- 
haloes have 0.65 < 02/01,03/02 < 0.75. Those shown in 
the lower left panel have 0.55 < 02/01,03/02 < 0.65, and 
those in the lower-right panel are highly fiattened systems, 
with 0.45 < 02/01,03/02 < 0.55. For comparison, we also 
plot in each panel the collapse barriers associated with the 
Bond et al. (1991) ellipsoidal collapse model (dot-dashed 
line), the Sheth et al. (2001) fit (dashed line), and the "true" 
ellipsoidal collapse barrier consistent with the quoted proto- 
halo shape (solid curve) . As in Figure 7, we mark each curve 
with the percentage of haloes that lie above it. 

There are a couple of things in Figure 8 worth high- 
lighting. First, note that the distribution in {eh,5h) at fixed 
halo m,ass depend sensitively on the the intrinsic shape of 
the proto-halo. For example, of objects with 02/01 and 
03/01 > 0.8 all have e ^ 0.42. For the strongly triaxial 
proto-haloes this distribution is shifted to the right; for ex- 
ample, we find that only 21.8% of all proto-haloes with 
0.45 < 02/01,03/02 < 0.55 have e < 0.42. For all haloes 
in this mass bin only about half have e < 0.42. The reason 
for this is clear in the context of ellipsoidal collapse dynam- 
ics. Perturbations that are initially very elongated require 
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Figure 8. Same as Figure 7 but for proto-haloes in the mass range 24.3 < MFOF/(10^''/i^^AfQ) < 48.7. Each panel plots the subsample 
of haloes in a narrow bin of 02 /ai ~ a:i/a2 (indicated in each panel). Dot-dashed and dashed lines are the same as those plotted in 
Figure 7; solid curves show the collapse threshold obtained from the ellipsoidal collapse model assuming an initially triaxial perturbation 
with 02/01 = 03/02 chosen to match the haloes in each subsample. 



external tidal compression along their major axis to balance 
the suppressive internal shear. Increasing the ellipticity of 
the tidal field in our model results in larger external com- 
pression along the long axis of the ellipsoid and ultimately 
reduces the density contrast required for collapse. 

Another important point to note is that, for any given 
proto-halo shape, the minimum barrier heights are described 
surprisingly well by our simple modification to the classic 
ellipsoidal collapse model of Bond et al. (1991). The nearly 
spherical proto-haloes plotted in the upper left panel of Fig- 
ure 8 are essentially always above the barrier inferred for 
the collapse of a spherical region in the presence of an exter- 
nal tidal field; indeed, the maximum ellipticity allowed for 
spherical proto-haloes is ~ 0.53 (for p — 0), for e greater 
than this value collapse along the minor axis of the tidal 
field is impossible, and so "full collapse" never occurs. Re- 
laxing the assumption that the proto-halo should be spheri- 
cal has a substantial impact on the inferred collapse bound- 
ary. For example, of the proto-haloes in this mass range with 
0.65 < 02/01, 03/02 < 0.75 nearly 44% have average internal 
density contrasts below what is required by the Bond et al. 
(1991) model for their measured ellipticities. Indeed, for the 



proto-haloes with the most extreme axis ratios (lower-right 
panel) we find that only about one in five come from regions 
of the linear density field that meet the collapse requirement 
based on the "standard" ellipsoidal collapse model, starting 
from initially spherical Lagrangian patches. 

These results will have important implications for an- 
alytic predictions of the halo mass function and clustering 
properties based on the excursion set method, and may help 
improve our understanding of structure formation in gen- 
eral. However, in order to be useful one must have a way to 
anticipate the typical shape of a proto-halo that is expected 
to form at a given position, smoothing scale, etc. We discuss 
this next. 



6 DISCUSSION 

In the ellipsoidal model the collapse threshold depends ex- 
plicitly on the shape of the tidal field at a given point x; 
the inferred collapse threshold can then be compared with 
the true overdensity at that location to decide whether or 
not collapse should take place. The approach of Sheth et al. 
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Figure 9. The fraction of haloes of a given mass whose smoothed 
linearly extrapolated overdensities exceed the threshold for col- 
lapse. Different colored curves show results for our two simula- 
tions; red lines for our 1200 h~^ Mpc box simulation; blue curves 
to our 150 h~^ Mpc box run. Dashed lines (with triangles) adopt 
the mass-dependence of the standard ellipsoidal collapse barrier 
of Sheth & Tormen; solid curves (with circles) adopt a collapse 
barrier computed independently for each mass bin assuming an 
initial density perturbation with axis ratios consistent with the 
average shape of proto-haloes in each bin. The vertical lines mark 
the 500 particle limit, below which shape estimates may be biased 
by low particle numbers. 



(2001) was to adopt, for every spatial location, the most 
probabilistic values of e and p in Gaussian random fields, 
resulting in a barrier height that depends explicitly on filter- 
ing scale: Cmp = a{M)/5/\fE (since pmp = for any cr(M)). 
The resulting barrier height depends on only one variable: 
the characteristic mass scale under consideration. 

We have shown that relaxing the assumption that the 
assumption of sphericity for the collapsing Lagrangian re- 
gions leads to very different predictions for the collapse bar- 
riers at any e, p, or mass scale. But what is the characteristic 
shape of a proto-halo and how is it related to the exter- 
nal tides, {e,p), the characteristic mass scale, the filtering 
choice, etc? In this section we describe two simple methods 
for estimating in advance the potential shapes of collapsing 
proto-haloes at any position x. 

One possibility is to estimate the characteristic shape of 
collapsing proto-haloes directly from the mass dependence of 
proto-halo shapes shown in Figure 1. The best-fit relations 
given in eq. 7 can then be used to determine the average 
proto-halo shape expected for a given filter volume (which 
corresponds to a well-defined halo mass scale) . At any given 
smoothing scale, a{M), there is then a well-defined bound- 
ary 5b (e, p, M) that can be used to select possible sites for 
collapse based on the properties of the tidal field at x. 

In general, this requires a calculation of the full ellip- 



soidal collapse barrier at every characteristic mass scale, 
since the axis ratios are themselves smooth functions of 
mass. In order to make the calculation simpler we have as- 
sumed that p = for all x (so that the tidal field is fully char- 
acterized by its ellipticity, e) and that dark matter proto- 
haloes can be approximated as "perfectly tria:xial ellipsoids" , 
i.e. a2/ai = azjai^ so that a single axis ratio is sufficient for 
describing their shape. 

In Figure 9 we plot the fraction of proto-haloes whose 
linear density contrast exceeds the threshold for collapse 
based on the method outlined above. Here we estimate, at 
any given halo mass scale, the typical value of a2/a\ from 
eq. 7, set a^/ai — (a2/a{f' , and compute the corresponding 
barrier as a function of tidal ellipticity in the usual fashion, 
by assuming p = 0. At any given mass scale, the fraction 
of points that lie above this threshold are shown as solid 
circles connected by lines. Triangles connected by dashed 
lines show the fraction of haloes that lie above the barrier 
advocated by Sheth et al. (2001). We use blue and red lines 
to distinguish results obtained from our 150 h~^ Mpc and 
1200 h~^ Mpc box runs, respectively. 

Taking into account the average mass-dependence of 
proto-halo shapes we find that ^ 90% now lie above the 
ellipsoidal collapse boundary. Adopting the Sheth et al. 
(2001) barrier instead yields very different results: at ~ 
lO^^/i~^M0, for example, only about 70% of proto-haloes 
should collapse given their averaged ellipticities and density 
contrasts. At the highest mass scales, where haloes are ini- 
tially more spherical and have low ellipticities, roughly one 
half of proto-haloes should not collapse in the standard ap- 
proach to ellipsoidal collapse. This is a strong indication that 
carefully accounting for proto-halo shapes (and ultimately 
their boundaries) is key to realistic excursion set modeling. 

Another possibility is to attach potential proto-halo 
shapes to points in the smoothed linear density field by ex- 
ploiting the correlations between the ratios ai/aj and Xi/Xj 
of the eigenvalues of the inertia and tidal tensors (see Fig- 
ure 3). After some experimentation we found a relatively 
well-defined and mass-independent trend between the ratio 
as/ai and A3/A1; we show the median trends and their as- 
sociated one-sigma scatter for several different mass bins in 
Figure 10. Note that the haloes selected for this plot span 
more than four orders of magnitude in mass, yet at any given 
mass scale the average minor-to-major axis ratios, aa/ai, are 
highly correlated with the relative compression factors along 
the minor and major axes of the tidal field. The main trend 
can be reasonably well described by the simple function 



a-i 



C X tanh , 

Hi \ Al 



A3 



-V', 



(25) 



where C — 0.465 and tp = 0.488 are the best-fit parameters. 
Eq. 25 is shown as a dotted-dashed line in Fig. 10 

Given eq. 25, one can estimate the axis ratios as/ai of 
potential proto-haloes completely in terms of the tidal field 
measured at a given point, without making any explicit ref- 
erence to mass scale. Given an estimate of azjax we can set 
a^ijo,! = -v/aa/ai (i.e. assume perfectly triaxial proto-haloes) 
and compute the shape-dependent collapse boundaries on a 
point by point basis with knowledge of only the eigenvalues 
of the linear tidal field: Ai > A2 > A3. The collapse barrier 
can then be compared to the true overdensity at that point 
in the same sense as before. 



14 Ludlow & Porciani 



1,0 



0.6 - 



0.4 



0.2 



M/(10'° K' Ms) 

(24.3 - 48.7) 

(60.8 - 121.7) 

(243.3 - 486.6) 
(12461 - 24921 
■(31152 - 62303' 
(62303 - 124607) (1372) 



(20788 
9268) -mr 

(2786 
(24204) h- 

f5777 ^ 




-1.0 



-0.5 



0.0 



0.5 



1.0 



1.2 


1 1 1 


1 < < , 1 < , , 1 < , , 1 , , < 1 , , 






^- 


--ASheth, Mo & Tormen (2001) 






• 


— •Ellipsoidal Collapse 




1.0 




T T 


- 










y^ y^Zr^-^'^'^^ 


: 






•^^ / -^ 




0.8 


{ 


( ^f-^ 


- 




/ 


M ^ -^ 


- 


o 


- /. 


^ /^ ^'^ 


- 


t 0.6 

to 


4 


/ ^ 


- 






/ / 


. 






f X 


. 


0.4 


1 / 

■ / / 


- 




/ 


X 






/ 


/ 




0.2 


- i 


/ 


- 






i 


- 






/ 


- 






X 


. 


0.0 


1 1 1 


1 1 1 . 1 1 11 1 1 11 1 11 1 1 . . 





-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 



Figure 10. Median values of axis ratio aa/ai for proto-haloes 
computed in bins of A3/A1. Different colored curves and symbols 
denote different mass bins, as indicated in the legend. The total 
number of haloes in each bin is also given. The error bars cor- 
respond to the 25*'' and TS"* percentiles. There is a strong and 
relatively mass independent trend between 03 /ai, which can be 
approximated by eq. 25, as shown with the triple-dot-dashed line. 



In Figure 11 we plot the fraction of proto-haloes whose 
average overdensities exceed the collapse threshold inferred 
in this way. Haloes are binned according to the ratio A3/A1 
measured at their centers of mass and the appropriate col- 
lapse threshold is then computed from eq. 25. The curves 
shown in Figure 11 adopt the same colors and line-styles as 
those in Figure 9. Solid curves with circles correspond to the 
fraction of proto-haloes with measured density contrasts ex- 
ceeding the barrier height inferred from our model; dashed 
curves to the fraction of haloes with & > Sec, the Sheth 
et al. (2001) model (eq. 2). Clearly, more careful modeling 
of the intrinsic shapes of dark matter proto-halos can bet- 
ter account for the large scatter in their measured density 
contrasts. 



7 SUMMARY 

We have used two cosmological simulations of structure for- 
mation to test a key assumption commonly used in excursion 
set modeling of structure formation, namely, that the regions 
of the linear density field that collapse to form bound objects 
are spherical. We used our results to motivate a deeper look 
into the collapse boundaries of dark matter halos predicted 
by the ellipsoidal collapse model. Our main results can be 
summarized as follows. 

• The shapes of the Lagrangian patches that eventually 
collapse to form CDM halos are highly non-spherical. Over 
the mass scale ~ lO^/i'^M© to ~ lO^^/i'^M© the average 
intermediate-to-major axis ratio varies from {a2/ai) ~ 0.7 



Figure 11. The fraction of haloes with smoothed, linearly ex- 
trapolated overdensities above the threshold predicted by the el- 
lipsoidal collapse model. Fractions are computed in bins of A3/A1 , 
which measures the relative compression factors along the minor 
and major axes of the linear tidal field. Boundaries are computed 
from the ellipsoidal collapse model described in Section 4, with 
proto-halo shapes estimated from eq. 25. In the model, all proto- 
haloes are assumed to be perfectly triaxial, and we set p = for 
simplicity. The solid curves plot the fraction of haloes with mea- 
sured density contrasts exceeding the collapse threshold inferred 
from our model; dashed curves plot the fraction of haloes exceed- 
ing the model of Sheth et al. (2001). Different colored curves show 
results for our two simulations; red lines for our 1200 h~^ Mpc 
box simulation; blue for our 150 h~^ Mpc box run. 



to ~ 0.8, and {a^/ai) from roughly ^ 0.4 to ~ 0.65. 
Above the characteristic mass scale for collapse (roughly 
W^^h~^MQ) proto-haloes are, on average, "perfectly triax- 
ial", with axis ratios a^/ai ~ 03/02. 

• The shapes of dark matter proto-haloes are primarily 
determined by the relative compression factors associated 
with the linear tidal field acting upon them (see Fig 3). 
For example, haloes that form in regions of the linear den- 
sity field where Ai is large and positive, but A3 is large and 
negative, tend to accrete material from very fiattened and 
elongated regions. This is because material is compressed 
along one axis of the tidal field, and dilated along an or- 
thogonal direction. The net effect of this push-and-puU of 
material is a strong alignment between the principle axes of 
the inertia tensors of proto-haloes and those of the tidal field 
measured at their geometric centers (Porciani et al., 2002; 
Lee & Pen, 2000). The origin of this correlation is rooted in 
the tidal field, which sets the boundaries of regions of the 
linear density field that are able to collapse by a given time. 

• The alignment between the eigenvectors of the linear 
tidal field and the inertia tensor of proto-haloes can be ex- 
ploited to simplify the application of ellipsoidal collapse to 
triaxial density perturbations. In Section 4 we modified the 
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standard ellipsoidal collapse model of Bond et al. (1991) 
in order to account for the non-spherical nature of the La- 
grangian patches that collapse to form halos. Relaxing the 
assumption that these regions are initially spherical yields 
very different predictions for the thresholds for non-linear 
structure formation. 

• In the absence of external tidal forces any non-spherical 
density perturbation requires a higher density contrast for 
collapse than the corresponding spherical perturbation. This 
is because the internal tidal forces suppress collapse along 
the major axis of the perturbation, and correspondingly 
higher densities are needed to counter this. When external 
tidal forces are introduced as described in Section 4, there is 
a well-defined minimum barrier height for any non-spherical 
perturbation; for a given tidal prolateness, p, this minimum 
occurs when the density contrast for collapse along the ma- 
jor and minor axis of the density perturbation are equal (See 
Fig 5). 

• We provided two possible methods for estimating the 
typical shape of a proto-halo that may form at any given lo- 
cation in the smoothed linear density field. The first involves 
measuring the average shapes of proto-haloes as a function 
of their mass, and using these shapes to compute explic- 
itly the average barrier height for gravitational collapse as 
a function of halo mass. Another possibility is to estimate 
directly the proto-halo shape by exploiting the correlation 
between the axis ratio as/ai and the relative tidal com- 
pression factors A3/A1 (see Fig. 10). Since proto-halos are 
roughly triaxial (03/02 = a^/ai) one can compute the typi- 
cal shape of a proto-halo that will form at a given location 
based solely on the ratio A3/A1 measured at that location. 
As many as ^ 90% of (A''fof > 500) dark matter proto- 
haloes have measured linear density contrasts that exceed 
the collapse thresholds inferred in these ways. 

All in all, our results indicate that the standard ellip- 
soidal collapse model of Bond et al. (1991) provides a rather 
incomplete census of the possible collapse thresholds for the 
formation of CDM halos from Gaussian random fields. In 
addition to the specifics of the linear tidal field at a given lo- 
cation, the collapse threshold depends on the assumed shape 
of the Lagrangian region that eventually collapses to form 
a bound object. This is because the linear overdensity re- 
quired for collapse in a given tidal field is sensitive to the 
initial displacement of a given mass element from the per- 
turbation center, thus changing how the collapse barriers de- 
pend on auxiliary properties of the tidal field. These results 
will have important implications for excursion set modeling 
of the halo mass function, spatial clustering and mass accre- 
tion histories, and may provide new insight into the complex 
process of dark matter halo formation. 
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